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O Abstract 

^ Using simulations of hard rods in smectic-A states, we find non-gaussian diffusion and het- 

S erogeneous dynamics due to the equilibrium periodic smectic density profiles, which give rise to 
permanent barriers for layer-to-layer diffusion. This relaxation behavior is surprisingly similar to 

O 

O that of non- equilibrium supercooled liquids, although there the particles are trapped in transient 

^ (instead of permanent) cages. Interestingly, we also find stringlike clusters of up to 10 inter-layer 

^ rods exhibiting dynamic cooperativity in this equilibrium state. 
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Smectic liquid crystals consist of stacks of fluidlike layers of oriented rodlike particles [T] . 
They can be stabilized by attractive Van der Waals forces, e.g. in thermotropic systems of 
mesogenic molecules [T] , or by repulsive hard-core interactions between sufficiently elongated 
particles [2j . In the relatively simple smectic-A phase the static structure is characterized by 
long-range orientational ordering of the rods combined with a one-dimensional periodic den- 
sity variation in the direction parallel to the rods. While the equilibrium properties of smectic 
phases are relatively well understood [1], little is known about their dynamics on the particle 
scale, even though Helfrich's early report of diffusion ("permeation") of anisotropic particles 
through smectic layers goes back 40 years [3]. Recently, however, exciting progress was made 
based on newly developed experimental techniques (e.g. NMR coupled to strong magnetic 
field gradients [1], or fiuorescent labelling of rods [5j), which revealed direct observations 
of non-Gaussian diffusion and quasiquantized layer-to-layer hopping across a barrier [5]. 
This triggered new theoretical work, based on dynamic density functional theory (DDFT), 
which not only confirmed the non-Gaussian diffusive motion and the one- dimensional "per- 
manent" barriers due to the static smectic background, but also showed the importance of 
"temporary" cages formed by neighboring rods [6]. 

Interestingly, non-Gaussian diffusive behavior due to heterogeneous dynamics of "slow" 
and "fast" particles is also a key feature in glassy systems and supercooled hquids, in which 
individual particles are trapped in transient cages formed by their neighbors. This het- 
erogeneous dynamics explains the experimentally observed nonexponential relaxation and 
stretching of time correlation functions and is found to be closely related to cage rearrange- 
ments and to cooperative motion, in which a small fraction of the particles (typically a few 
percent) move collectively in stringlike [3 |8] or compact [9l [10] clusters . Several models 
based on a heterogeneous distribution of diffusion coefficients, jump times or jump lengths 
have been introduced to explain the heterogeneous dynamics due to temporary cages [TT] . 
while another model explains the non-Gaussian diffusion by the dynamics of a single Brown- 
ian particle in a periodic external potential, which is very similar to the permanent barriers 
of the smectic-A phase fl2]. The intriguing question that we address in this Letter is to 
what extent the dynamics in the smectic-A phase is heterogeneous and/or collective, or, in 
other words, to what extent is the equilibrium smectic-A dynamics resemblant to that of 
out- of- equilibrium quenched supercooled liquids? This issue cannot be addressed directly 
in fluorescence experiments, in which only a small fraction of the rods is labelled such that 
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moving clusters cannot be observed. While DDFT yields average quantities such as density 
profiles, Van Hove functions, and single-particle barriers, but does not provide cluster infor- 
mation [B]. We therefore resort to computer simulations of hard-rod fluids in the smectic 
phase. Contrary to the essentially one-particle analysis of Ref. [13], we do find clear evidence 
for heterogeneous dynamics and cooperative permeation of strings of 1-10 rods. 

We perform kinetic Monte Carlo (KMC) simulations in the isobaric-isothermal ensemble 
(NPT) of N = 1530 — 3000 freely rotating hard spherocylinders with aspect ratio L* = 
L/D = 5, where D is the diameter and L + D the length of the rods. For L* = 5 the smectic 
bulk phase melts into a nematic phase below P* = 1.4 and freezes into a crystal above 
P* = 2.3, where P* = PD^/keT is the reduced pressure with ks Boltzmann's constant 
[E]. We study the dynamics at two state points in the bulk smectic phase characterized by 
Pi = 1.6 and — 2-0, which correspond to packing fractions rji = 0.508 and 7]2 = 0.557, 
respectively. Translational and rotational moves are attempted and accepted if no overlap 
is detected. A rectangular simulation box with 5-10 smectic layers and periodic boundary 
conditions are employed. Volume changes are attempted every MC steps by randomly 
changing the side length of the simulation box. In KMC simulations it is convenient to 
relate the number of MC cycles to time t. We chose r = D^/Dtr as our unit of time, where 
Dtr is the translational short-time diffusion coefficient, which is the isotropic average of the 
diffusion coefficients in the three space dimensions. 

Denoting the positions of the rods by rj = (xj, y,, Zi), we first measure the (relative) prob- 
ability 7r(z) of finding a rod at position z, where the 2;-axis is parallel to the nematic director 
n. Following Ref. [5] we introduce the Boltzmann factor 7t{z) oc exp (— f/ {z)/kBT) with U{z) 
the effective potential for diffusion out of the middle of a smectic layer. Fig. 1 shows U (z) 
in a small z-regime for our two state points together with the fit U{z) = Uq sin^ (irz/h), with 
barrier height Uq = (3.5, 7.5)kBT and layer spacing h = {L + D){1.05, 1.03) as fit parameters 
for the two states. The denser state thus reveals a substantially larger diffusion barrier and 
a stronger confinement to the middle of the smectic layers. 

In order to further quantify the physical picture of a hopping-type ^-diffusion, with the 
rods rattling around in a given layer until they overcome the free-energy barrier and jump 
to a neighboring layer, we calculate the self part of the Van Hove correlation function 
(VHF) Gs{z,t) = ^ (^J2iLi {zi{tQ + t) — 2j(to))]^, with (...) the ensemble average over 
all particles and initial time to, and 6 the Dirac-delta. Note that Gs{z, t) gives the distribution 
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FIG. 1: Effective potential U {z) at two state points in tlie bulk smectic phase with packing fractions 
r/i = 0.508 (red squares) and rj2 = 0.557 (black circles); the lines are fits. The horizontal bars denote 
the standard deviation a = (0.6,0.3)D at r/i and 772, respectively. 

for the particle 2;-displacements during a time-interval t, and would be a Gaussian of z for 
freely diffusive particles. Fig. 2(a,b) show the VHF for our two state points as a function 
of z for several equidistant t's, showing the appearance of peaks at integer layer spacings 
consistent with earlier experimental [S] and theoretical [B] results. The height and spatial 
extension of the peaks is larger at r/i, indicative of a faster layer-to-layer diffusion than at 772- 
Non-gaussian VHF's were analyzed in 2D liquids [15], glasses [IT], and liquid crystals [HIE], 
and attest to the presence of dynamical heterogeneities which imply particles diffusing faster 
or slower than the average. This is quantified by the well-known non-gaussian parameter 
(NGP) [IS], for z-diffusion defined by a,(t) = -1 + {Az^tf) m^z^iff with Az{t) = 
z{tQ + t)—z(tQ) the 2;-displacement of a rod in the time interval t starting at to- Likewise non- 
Gaussian in-plane diffusion can be characterised by an NGP axy{t). Heterogeneous dynamics 
occurs on a time-scale t if az{t) or axyit) is non-vanishing. In Fig. 2(c,d) we show the NGP's 
for our two state-points, with az{t) showing a clear peak at t* ~ (2, 10)r, while <yxy{t) is 
hardly visible on the scale of the figure and thus essentially vanishes for all t. In other words. 
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FIG. 2: (a,b) Self part of the Van Hove function Gs{z^ t), (c,d) non-Gaussian parameters a^(i) and 
(^xyit), (e,f) mean-square displacement (MSD) in units of D^, and (g,h) self-intermediate scattering 
functions Fs{t), at packing fraction rji = 0.508 (left column) and 772 = 0.557 (right column). The 
time intervals in (a,b) range from t = OAt (dotted lines) to t = 40r (solid lines), with increments 
of 8r. The solid and open symbols refer to the z and xy direction, respectively. 



in-plane diffusion is dynamically homogeneous while permeation is heterogeneous, with t* 
and az{t*) increasing with density. In order to determine the relaxation dynamics we plot 
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the mean square displacements (MSDs) (A^^^^)^ ^nd (Ax2(t) + A?/2(t)) in Fig. 2(e,f) for the 
two state points. The xy-MSD shows a smooth crossover from short- to long time diffusion, 
while the 2;-MSD develops an intermediate cage-trapping plateau up to t* beyond which 
interlayer diffusion takes place. Note that the 2;-MSD exceeds the xy-MSD only at short 
times. 

In addition we study structural relaxation by computing the self-intermediate scattering 
function Fs{t) = (exp[zq ■ Ar(t)]) at wavevectors qD = {0,0, q^) and {qx,(ly,0), with ^2 = 1 
and (g^ + Q^Y^'^ = 6, that correspond to the main peaks in the static structure factor. Here 
Ar(t) denotes the particle displacement during a time-interval t. Results are shown in Fig. 
2(g,h). The in-layer dynamics is very fast with nearly exponential decay, typical of fluid-like 
behavior. By contrast, the inter-layer dynamics is much slower, with Fs{t) decaying in two 
steps at t' and t" separated by a plateau during t' < t < t" . Interestingly this plateau 
(corresponding to /3-relaxation) coincides with the nonvanishing of az{t), and can hence be 
attributed to the heterogeneity of the dynamics [7]. The first step at t'/r < 1 corresponds 
to the rattling of rods inside the temporary cage and permanent smectic background formed 
by neighboring rods [6j, whereas the second one (a-relaxation) corresponds to the escape 
on a time scale that increases from t" /t ~ 100 at rji to 10^ at ri2. The increase of t" and 
that of the height of the plateau of Fgif) with density is also observed for colloidal glasses 
[T7] . At rji, the long-time decay of Fs{t) is well fitted by a stretched exponential of the 
form exp[— (t/t^))'^], with relaxation time t^/r = 650 and (3 = 0.6, once more confirming 
the heterogeneous nature of the relaxation dynamics |T7j; tr at 772 is beyond our simulation 
time. 

Having established the heterogeneity of the dynamics in the equilibrium smectic bulk 
phase, our next question concerns the identification of the fast particles. Comparing the 
VHF at t* of Fig. 1(a) and (b) with a gaussian approximation with the same MSD reveals 
that the fast-moving particles, which are responsible for the non-gaussian character, have 
travelled over more than h/2 during a time interval t*. Therefore they are intimately related 
to interlayer particles, which reside more than some distance 6 from the nearest smectic 
plane. In order to define 6 sensibly, we consider the variance of one period of 7t{z) as 
^2 ^ j'^l^^^z'^n{z)dz, which gives a/D = (0.6,0.3) for rji and 772, respectively, as indicated 
in Fig.l. We now set 6 = ka with k = 1,2,3. The fraction of so-called interlayer 
particles is then /i = (0.28,0.45), /a = (0.05,0.13) and /a = (0.02,0.03) for our two state- 
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points, showing, perhaps surprisingly, that the denser state contains more interlayer particles 
caused by the smaller o. We also calculated the pair distribution of the interlayer particles 
(not shown), revealing a higher contact value than that of the bulk smectic phase, suggesting 
substantial clustering of interlayer particles despite their low concentration /fc?]. Stringlike 
clusters composed oi n = 1, ■ ■ ■ ,10 interlayer rods can indeed be identified, as shown for 
5 = 2cr by the size distribution P{n) in Fig. 3. Our (rather stringent) cluster criterion is 
such that two interlayer particles belong to the same cluster if the z and xy distance are 
smaller than h and D, respectively. At rji, one deduces from P{n) that ~ 95% of the clusters 
consist mostly of 2 or 3 rods, while clusters of more than 5 rods are rare but do exist. The 
denser smectic phase has larger clusters, which is again similar to supercooled liquids and 
glassy systems, in which the cluster size increases with increasing cage-trapping [HI Hn] • The 
fit P{n) oc exp{—an) is accurate with a = (1.5,0.7) at rji and 772, respectively, from which 
the average cluster size follows as {n) = (1 — exp(— a))^^ These results depend, however, 
strongly on 6 as revealed by the inset of Fig. 3. 

With all these resemblances to glassy dynamics we now study the dynamics of the layer- 
to-layer diffusion, for which Fig. 4b shows some typical trajectories. Most of the rods 
"jump" fast compared to the dwelling time within a smectic layer. Some rods diffuse to 
the inter-layer spacing and return to their original layer (yellow trajectory), others move 
from one layer to another several times (blue), in some cases the dwelling time in the inter- 
layer spacing is quite long (cyan), while double jumps can be also observed (pink). These 
observations suggest a rather broad distribution of jump-times tj to diffuse from one layer 
to the next, where tj is defined as the time span between first and last "contact" with the 
new and old layer, respectively, with contact established if the rod is at a distance S from 
the middle of a smectic layer. At rji and for 6 = 2cr, the distribution of jump-times n(tj) 
as obtained by averaging over many trajectories is plotted in Fig. 4a, and shows that the 
most-probable jump-time is 0.14r, the median is at t*j = 0.2r, while the distribution extends 
over two decades 0.01 < tj/r < 1. Having developed a sense of time, we can characterize the 
trajectories further. We distinguish single and multiple jumps, depending on the dwelling 
time in the next layer being longer or shorter, respectively, than t*j before another jump 
is started. At 771 we find that ~ 5% of the total number of jumps is multiple, of which 
~ 83% is a double and the remaining fraction a triple jump. At 772? where t*j = 0.27r, 
the fraction of multiple jumps is 0.6%. The inter-layer rods are generally oriented along n. 
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FIG. 3: Probability distribution P(n) of the number n of inter-layer rods in a stringlike cluster 
(with 6 = 2a, see text), at r/i (#) and ri2 (O). The inset shows P(n) at r/2 for 6 = a (□), 2a 
(O), and 3a (O). The solid lines denote the fit P(n) oc exp(— an) with a given in the figure. The 
snapshot of 3000 rods at r]i shows predominantly in-layer rods (green) and single interlayer rods 
(black), both with diameters reduced to D/4 for clarity. The thicker rods denote transverse ones 
(red), as well as clusters of 2 (yellow), 3 (orange), 4 (pink), 5 (blue), and > 6 (brown) rods, (color 
online) . 



Remarkably, however, at rji also transversely oriented rods in between two smectic layers 
have occasionally been observed [HI [19], which diffuse either to a new layer or return to the 
original one by rotating parallel or anti-parallel to their original orientation. 

With our clear evidence for heterogeneous dynamics and spatial correlations between 
the fast-moving interlayer particles, it seems natural to investigate dynamic cooperativity, 
i.e., whether or not the stringlike clusters observed in static configurations actually move 
collectively. This requires a cluster criterion that not only involves spatial but also temporal 
proximity to identify collectively moving rods. Two jumping rods i and j are considered to 
move cooperatively if their arrival times t'-'^ and t^^^ in their new layers (i.e. the first time 
at which their distance to the middle of the new layer equals 5) satisfies \t^'^^ — t'^-'-'l < At, 
while rj(t'^*)) and Vj{t^^^) satisfy our static spatial cluster criterion. Using 5 = 2a and 
At = r (i.e. long enough for any jump to finish according to Ii{tj)), we find at rji that 
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the number of collective jumps is ~ 25% of the total number of jumps, and involves mainly 
2 79%) or 3 17%), and very rarely > 4 (< 4%) rods. Interestingly, 42% of the 
collective jumps involve rods diffusing in opposite directions. These characteristic values are 
insensitive to small modifications of our spatial cluster criterion, while a smaller temporal 
interval of At = t*j (the median jump time) reduces the fraction of collective jumps to 
~ 9%. In other words, the motion is indeed strongly cooperative at r]i. If the analogy 
with glassy systems would hold even further, one would expect more cooperativity at ■ 
Surprisingly, perhaps, despite the larger static clusters we find less collective motion at 
r]2. For instance, with 6 = 2(7 and At = 1.2r (maximal jump time) only ^ 4% of the 
jumps can now be regarded as collective (involving essentially only 2 rods). We attribute 
this reduction of cooperativity upon approaching crystallization to the higher permanent 
barriers that reduce the probability of static clusters to actually complete their attempted 
jumps. Such permanent barriers do not exist in undercooled liquids (where the dynamic 
cluster size grows upon approaching the glass transition |H1 [ID]), whereas they do exist in 
crystals, where structural relaxation is virtually non-existent. 




tj/T x/D 

FIG. 4: (a) Distribution n(t j) of layer-to-layer jump times tj, based on S = 2a at density r]i. (b) 
Trajectories of jumping rods projected onto the xz plane, with the dashed lines representing the 
middle of the smectic layers (color online). 

In conclusion, for equilibrium hard rods in the smectic phase we find non-gaussian and 
cooperative relaxation behavior that is surprisingly similar to that of non- equilibrium super- 
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cooled liquids. We attribute the non-gaussian diffusion to the spatial inhomogeneous density 
of the smectic phase, and connect this to djTiamic heterogeneities and to static clusters of 
inter-layer rods. Interestingly, we also find clear evidence for cooperative motion of stringlike 
clusters. Our results might be relevant for dynamics in other inhomogeneous liquids, e.g., 
confined fluids (micro and nanofluidic devices) and columnar liquid crystals. 

We thank P. van der Schoot and M. Bier for useful discussions. This work was financed 
by a NWO-VICI grant. 
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